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Abstract — An adaptive randomized distributed space-time cod- 
ing (DSTC) scheme is proposed for two-hop cooperative MIMO 
networks. Linear minimum mean square error (MMSE) receiver 
filters and randomized matrices subject to a power constraint 
are considered with an amplify-and-forward (AF) cooperation 
strategy. In the proposed DSTC scheme, a randomized matrix 
obtained by a feedback channel is employed to transform the 
space-time coded matrix at the relay node. The effect of the 
limited feedback and feedback errors are considered. Linear 
MMSE expressions are devised to compute the parameters of 
the adaptive randomized matrix and the linear receive filters. 
A stochastic gradient algorithm is also developed with reduced 
computational complexity. The simulation results show that the 
proposed algorithms obtain significant performance gains as 
compared to existing DSTC schemes. 

I. Introduction 

Cooperative multiple-input and multiple-output (MIMO) 
systems, which employ multiple relay nodes with antennas 
between the source node and the destination node as a dis- 
tributed antenna array, apply distributed diversity gain and 
provide copies of the transmitted signals to improve the 
reliability of wireless communication systems IfTZl . Among 
the links between the relay nodes and the destination node, 
cooperation strategies, such as Amplify-and-Forward (AF), 
Decode-and-Forward (DF), and Compress-and-Forward (CF) 
1 1 1 and various distributed space-time coding (DSTC) schemes 
in (21, (3) and |fl9ll can be employed. 

The utilization of a distributed STC (DSTC) at the relay 
node in a cooperative network, providing more copies of the 
desired symbols at the destination node, can offer the system 
diversity gains and coding gains to combat the interference. 
The recent focus on the DSTC technique lies in the design of 
delay-tolerant codes and full-diversity schemes with minimum 
outage probability. An opportunistic DSTC scheme with the 
minimum outage probability is designed for a DF cooperative 
network and compared with the fixed DSTC schemes in |@). 
An adaptive distributed-Alamouti (D-Alamouti) STBC design 
is proposed in for non-regenerative dual-hop wireless 
systems which achieves the minimum outage probability. 

The channel state information (CSI) is very important for 
a wireless communication system and can be estimated by 
sending a block of training symbols to the destination node. 
The feedback technique allows the destination node to transmit 
the CSI or other information back to the source node, in order 
to achieve gains by pre-processing the symbols. In [6], the 
trade-off between the length of the feedback symbols, which 



is related to the capacity loss, and the transmission rate is 
discussed, and in Q, one solution for this trade-off problem 
is derived. The use of limited feedback for STC encoding 
has been widely discussed in the literature. In 0, the phase 
information is sent back for STC encoding in order to maintain 
the full diversity, and the phase feedback is employed in (9) to 
improve the performance of the Alamouti STBC. The limited 
feedback is used in ifTUl and ifTTl to provide the channel 
information for the pre-coding of an OSTBC scheme. 

In this paper, we propose an adaptive linear receiver de- 
sign algorithm with randomized distributed space-time coding 
optimization based on the MSE criterion for cooperative 
MIMO relaying systems with limited feedback. We focus 
on how the randomized matrix affects the DSTC during 
the encoding and how to optimize the linear receive filter 
with the randomized matrix iteratively. It is shown that the 
utilization of a randomized matrix benefits the performance 
of the system compared to using traditional STC schemes. 
Then an adaptive optimization algorithm is derived based on 
the MSE criterion subject to constraints on the transmitted 
power at the relays, with the aid of a stochastic gradient (SG) 
algorithm in order to release the destination node from the 
high computing complexity of the optimization process. The 
updated randomized matrix is transmitted to the relay node 
through a feedback channel with errors, and the influence of 
the imperfect feedback is discussed. 

The paper is organized as follows. Section II introduces 
a two-hop cooperative MIMO system with multiple relays 
applying the AF strategy and the randomized DSTC scheme. 
In Section III the proposed optimization algorithm for the ran- 
domized matrix is derived, and the results of the simulations 
are given in Section IV and Section V leads to the conclusion. 

II. Cooperative System Model 

The communication system under consideration, shown 
in Fig.l, is a cooperative communication system employing 
multi-antenna relay nodes transmitting through a MIMO chan- 
nel from the source node to the destination node with feedback 
channels to the relay nodes. The 4-QAM modulation scheme 
is used in our system to generate the transmitted symbol vector 
s[i] at the source node. There are n r relay nodes with N anten- 
nas for transmitting and receiving, applying an AF cooperative 
strategy as well as a DSTC scheme, between the source node 
and the destination node. A two-hop communication system 
that broadcasts symbols from the source to n r relay nodes as 



well as to the destination node in the first phase, followed by 
transmitting the amplified and re-encoded symbols from each 
relay node to the destination node in the next phase. After 
decoding at the destination node, the information matrix for 
encoding will be quantized first, and then transmitted back to 
each relay node through a feedback channel with noise and 
interference. The relay nodes quantize the feedback symbols 
and use them as a part of the encoding matrix in the next 
transmission. We consider only one user at the source node in 
our system that has N Spatial Multiplexing (SM)-organized 
data symbols contained in each packet. The received symbols 
at the k — th relay node and the destination node are denoted 
as rsR k and r$D, respectively, where k — 1,2, ...,n r . The 
received symbols rsR h will be amplified before mapped into 
an STC matrix. We assume that the synchronization at each 
node is perfect. The received symbols at the destination node 
and each relay node can be described as follows 

r S R k [i] = F k [i] s [i] + n SR k [i] , ( 1 ) 

r SD [i] =H[i]s [i] + n SD [i] , (2) 

i = 1,2, ... ,N, k = 1,2, ... n r , 

where the N x 1 vector nsR k [i] and nsn[i] denote the zero 
mean complex circular symmetric additive white Gaussian 
noise (AWGN) vector generated at each relay and the destina- 
tion node with variance a 2 . The transmitted symbol vector s[i] 
contains N parameters, s[i) — [s±[i], S2[i], sjv[i]], which 
has a covariance matrix _E[s[i]s H [/]] = all, where E[] 
stands for expected value, (•) denotes the Hermitian operator, 
a 2 is the signal power which we assume to be equal to 1 and 
I is the identity matrix. Fk [i] and H[i] are the NxN channel 
gain matrices between the source node and the k — th relay 
node, and between the source node and the destination node, 
respectively. 

After processing and amplifying the received vector rs Rk [i] 
at the k — th relay node, the signal vector ssR k [i] = 
AR k D[i](Fk[i]s[i] + nsR k [i]) can be obtained and will be 
forwarded to the destination node. The amplified symbols 
in ssR k [i] will be re-encoded by an N x T DSTC scheme 
M(s[i]) and then multiplied by an N x N randomized matrix 
in [20], then forwarded to the destination node. The 
relationship between the k — th relay and the destination node 
can be described as 

R RhD \i] = G k \i]V\\i]MR kD \i] + N RkD \i], (3) 

k — 1 j 2. . ■ ■ . . 77ij- . 

where the W x T matrix M RkD [i] is the DSTC matrix 
employed at the relay nodes whose elements are the amplified 
symbols in ss Rk [i}. The N x T received symbol matrix 
RR k D[i] in (01 can be written as an NT x 1 vector rR k r>[i] 
given by 

iV 

r R k D [i] = \ i } G eg kj \i}ssR kj W + n RkD [i], (4) 

i=i 

where the NT x N matrix G eqkj [i] stands for the equivalent 
channel matrix which is the DSTC scheme M"(s[i]) combined 
with the channel matrix GR k D[i] and the block diagonal 
NT x NT matrix 0i eqkj [i] denotes the equivalent randomized 



matrix assigned for the j — th forwarded symbol at the 
relay node. The NT x 1 equivalent noise vector riR k D[i] 
generated at the destination node contains the noise parameters 
in NR k £)[i]. After rewriting RR k r>[i] we can consider the 
received symbol vector at the destination node as a N(n r + 1) 
vector with two parts, one is from the source node and another 
one is the superposition of the received vectors from each relay 
node, therefore the received symbol vector for the cooperative 
MIMO network we considered can be written as 

r[i]= Ef=l H e q M S M 

[ Y2UT,f=i^eq kj [i\G eqkj [i]s SRkj [i] _ 



N 

where the (T + 1)JV x (n r + 1)N block diagonal matrix 
D Dj [i] denotes the channel gain matrix of all the links in the 
network for the j — th symbol in [i] which contains the 
NxN channel coefficients matrix H[i] between the source 
node and the destination node, the NT x N equivalent channel 
matrix G eqk [i] for k — 1,2, ...,n r between each relay node 
and the destination node. The (n r + 1)N x 1 noise vector 
nfl[i] contains the received noise vector at the destination 
node and the amplified noise vectors from each relay node, 
which can be derived as an AWGN with zero mean and 
covariance matrix a 2 (l+ || 0\ eq k [i]G e qk [i]A RkD [i] \\ 2 F )I, 

where || X \\ F = \J^{X H ■ X) = Jtt{X ■ X H ) stands 
for the Frobenius norm. 

III. Joint Constrained Adaptive Randomized STC 
Optimization and Linear MMSE Receiver Design 

As derived in the previous section, the DSTC scheme used 
at the relay node will be multiplied by a randomized matrix 
subject to a power constraint before being forwarded to the 
destination node. In this section, we design a constrained 
adaptive optimization algorithm based on an SG estimation 
algorithm 1 1 8 1 for determining the optimal randomized matrix 
and the linear MMSE receive filters. 

A. Linear MMSE Receiver Design with RSTC Optimization 

The linear MMSE receiver design and the optimal RSTC 
matrices subject to a transmit power constraint at the relays 
are derived as follows. Define the (T + 1)N x 1 parameter 
vector Wj [i] to determine the j — th symbol Sj [i] . From (0 we 
propose the MSE based optimization with a power constraint 
at the destination node as 

[«>i[i].9W**]]=Mg min E[\\ Sj [i]-w?[i\r[i]f]., 
subject to 

N 

^trace(*H eqk .[i]< qk .[i])<P R , 
j'=i 

where r[i] denotes the received symbol vector at the desti- 
nation node which contains the randomized STC matrix with 
the power constraint of P R . If we only consider the received 



symbols from the relay node, the received symbol vector at 
the destination node can be derived as 

n r N 

k=l j=l 
n r N 

= E E ^ ® G <*., w^- w F i W a i w 
fe=ii=i 

n r AT ( 6 ) 

+ E E w Ge <? fcJ w a j w^si?., w 

+ n RD [i] 

n r N 

= E E ^ W s j w + W' 

fc=i J=i 

where C/.. [i] is an iVT x N matrix that contains all the 
complex channel gains and the amplified matrix assigned to 
the received symbol Sj[i] at the relay node, and the noise 
vector n u [i] is a Gaussian noise vector with zero mean 
and variance a 2 (l + Y%=i II ^■e qk] \i}G eqk] [i\A j \i] |||). 
Therefore, we can rewrite the MSE cost function as in (|7). 

Since Wj[i] can be optimized by expanding the righthand 
side of ^} and taking the gradient with respect to w* [i] and 
equating the terms to zero, we can obtain the j — th MMSE 
receive filter 



Wj \i] = {E[r\i}r H [i]}) E[r\i]sf\i]\ 



where E \r[i]r H [i]] denotes the auto-correlation matrix and 
E [r[i]Sj*[i]\ stands for the cross-correlation matrix. By op- 
timizing the randomized matrix 0i eqh [i] for each symbol 
at each relay node, we can first define a vector f eqk . = 
Ckj [i]sj[i] + Ckj [i]n.SRj, where the parameter nsRj denotes 
the j—th symbol in the noise vector nga, then the randomized 
matrix can be calculated by taking the gradient with respect 



to 9t egj( [i] and equating the term to zero, resulting in 



*W*] = KW(^ [f[z]f ff [z]] K-[z] + xi) 

\z,H 



(9) 



where E [f[i]r H [i]] denotes the auto-correlation matrix of the 
equivalent space-time coded received symbol vector without 
the randomized matrix at the relay node, and E \r H [i]sj[i]\ 
denotes the cross-correlation matrix. The power constraint can 
be achieved by multiplying the quotient of Pr and the trace 
of the updated randomized matrix. The expression in (|9j does 
not provide a closed-form solution of the randomized STC 
matrix 0i eqk . [i] assigned for the j — th received symbol at 
the k — th relay node because it requires the adjustment 
of the Lagrange multiplier A. This parameter needs to be 
adjusted in order to enforce the power constraint. Moreover, 
the expression in (|9]also requires an inversion calculation with 
a high computational complexity. With the increase of the 
number of antennas employed at each node or employing more 
complicated STC encoders at the relay nodes, the complexity 
increases exponentially according to the matrix size in (|9}- 



B. Adaptive Linear MMSE Receiver Design with Randomized 
Matrix Optimization Algorithm 

In order to reduce the computational complexity of the pro- 
posed design and compute the required parameters, an adaptive 
linear receiver design with randomized matrix optimization 
(ALRRMO) algorithm is proposed. We resort to a strategy 
that initially drops the power constraint, obtain the necessary 
recursions and then enforce the constraint with a normalization 
step. We define the Lagrangian of the constrained MSE 
minimization problem in (0 as 

Ciw^l^.m) = E -wf[i\r\i}\\ 2 ] 

N 



Etrace(?H cqkj [i]^ qkj [i])-P R J A, 

(10) 



J =1 



A simple adaptive algorithm for determining the linear receive 
filters and the randomized matrices can be achieved by taking 
the instantaneous gradient term of (0 with respect to w*[i] 
and with respect to 9i eqk respectively, which are 
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wf[i]r[i]) H r\ 
i]-wf\i]r[i\\ 



(11) 
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(8) where e 3 - [i] stands for the j — th detected error. After we obtain 



( fTTb the proposed algorithm is obtained by introducing a step 
size into the recursions. The proposed algorithm is given by 
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EjLi trace(5H oqkj [i + l]9lf qkj [i + 1] 



(12) 

where j3 and p, denote the step sizes for the recursions for 
the estimation of the linear MMSE receive filter and the 
randomized matrix in the RSTC scheme, respectively. The last 
equation in (fT2l) stands for the normalization of the randomized 
matrix after the iteration. According to (TTZt . the desired vector 
and the matrix depends on each other, so that the algorithm in 
ifTTl can be used to determine the linear MMSE receive filter 
and the randomized matrix iteratively, and the design can be 
achieved. The complexity for calculating the optimal Wj [i] and 
9\eq k] [i] is 0(N(T+1)) and 0(7V 2 T 2 ), respectively, which is 
much less than 0{2N 3 (T+ l) 3 ) and 0{2N 4 T 4 ) by using © 
and (0. As mentioned in Section I, the randomized matrix will 
be sent back to the relay nodes via a feedback channel which 
requires quantization as will be shown in the simulations. 

IV. Simulations 

The simulation results are shown here to assess the proposed 
scheme and algorithm. The system we considered is an AF 
cooperative MIMO system with the Alamouti STBC scheme 
using QPSK modulation in quasi-static block fading channel 
with AWGN, as derived in Section II. The bit error rate 



[WiM>9W,[i]] = ar § 



mm 



A- = l 3 = 1 



[i]C fc< [ 



i]sj[i] + n D [i] 



s.t. 



£ trace($H eqkj [i]9^ q [i]) < P F 



3=1 



(7) 



(BER) performance of the proposed adaptive linear receiver 
design with RSTC optimization algorithm is assessed, and the 
influence of the imperfect feedback channels are considered 
in the simulations. The system employs 1 relay node and 
each node in the system has 2 antennas. In the simulation, we 
define both the symbol power at the source node and the noise 
variance a 2 for each link to be equal to 1 . The RSTC scheme 
is designed by multiplying the 2x2 Alamouti STBC lfT5l by 
a randomized matrix with each element generated using e- 7 
where 9 is uniformly distributed in [0, 2w). 

The proposed ALRRMO algorithm is compared with the 
SM scheme and the traditional RSTC algorithm using the 
distributed- Alamouti (D- Alamouti) STBC scheme in [19| with 
n r = 1 relay nodes in Fig. 2. The number of antennas N = 2 
at each node and the effect of the direct link is considered. 
The results illustrate that without the direct link, by making 
use of the STC or the RSTC technique, a significant perfor- 
mance improvement can be achieved compared to the spatial 
multiplexing system. The RSTC algorithm outperforms the 
STC-AF system, while the ALRRMO algorithm can improve 
the performance by about 3dB as compared to the RSTC 
algorithm. With the consideration of the direct link, the results 
indicate that the cooperative diversity order can be increased, 
and using the ALRRMO algorithm achieves an improved 
performance with 2dB of gain as compared to employing the 
RSTC algorithm and 3dB of gain as compared to employing 
the traditional STC-AF algorithm. 

The simulation results shown in Fig. 3 illustrate the impact 
of the feedback channel for the ALRRMO algorithm. As 
mentioned in Section I, the optimal randomized matrix will 
be sent back to each relay node through a feedback channel. 
The quantization and feedback errors are not considered in 
the simulation results in Fig. 2, so the optimal randomized 
matrix is perfectly known at the relay node after the ALRRMO 
algorithm; while in Fig. 3, it indicates that the performance 
of the proposed algorithm will be affected by the accuracy of 
the feedback information. In the simulation, we use 4 bits to 
quantize the real part and the imaginary part of each element 
of the randomized matrix 9t e q k [£]> an d the feedback channel 
is a binary symmetric channel. As we can see from Fig. 3, 
by decreasing the error probabilities for the feedback channel 
with fixed quantization bits, the BER performance approaches 
the performance with the perfect feedback, and by making use 
of 4 quantization bits for the real and imaginary part of each 
parameter in the randomized matrix, the performance of the 
ALRRMO algorithm is about ldB worse with feedback error 
probability of 10~ 3 . 

In Fig. 4 and Fig. 5, the influence of different feedback 
error probabilities with various quantization bits are employed 
to test the performance of the ALRRMO algorithm. In Fig. 4, 
the BER performance with a perfect feedback channel is given 
as a lower bound with SNR = 15dB and SNR = 30dB, 
respectively. The error probability is fixed in P e = 10~ 3 . 
With the increase of the number of bits we employed in the 



quantization, the BER curves will approach the result with 
perfect feedback due to a more accurate estimation with the 
cost of computing complexity increase. With more quantiza- 
tion bits, the ALRRMO algorithm can achieve a performance 
as good as that with perfect feedback. However, it is worth to 
mention that the BER decreases slightly when we increase the 
number of quantization bits to 5 and 6. If we fix the number 
of quantization bits to 4, the BER performance gets worse 
with the increase of the feedback error probability as depicted 
in Fig. 5. By increasing the feedback error probability, the 
BER curves become gradually worse. Thus, there is a trade- 
off between the feedback error probability and the number of 
quantization bits. As indicated by the simulation results in Fig. 
3 to Fig. 5, the 4-bit quantization is an appropriate choice. 

V. Conclusion 

We have proposed an adaptive linear receiver filter design 
with randomized matrix optimization (ALRRMO) algorithm 
for the randomized DSTC in a cooperative system. A joint 
iterative estimation algorithm for computing the receive filters 
and the randomized matrix has been derived. The effect of the 
limited feedback and feedback errors are considered in the 
simulation. The simulation results illustrate the advantage of 
the proposed ALRRMO algorithm by comparing it with the 
cooperative network employing the traditional DSTC scheme 
and the fixed randomized STC scheme. The proposed algo- 
rithm can be used with different distributed STC schemes 
using the AF strategy and can also be extended to the DF 
cooperation protocol. 

References 

[1] J. N. Laneman and G. W. Wornell, "Cooperative Diversity in Wireless 
Networks: Efficient Protocols and Outage Behaviour", IEEE Trans. Inf. 
Theory, vol. 50, no. 12, pp. 3062-3080, Dec. 2004. 

[2] J. N. Laneman and G. W. Wornell, "Distributed Space-Time-Coded 
Protocols for Exploiting Cooperative Diversity in Wireless Networks", 
IEEE Trans. Inf. Theory, vol. 49, no. 10, pp. 2415-2425, Oct. 2003. 

[3] S. Yiu, R. Schober, L. Lampe, "Distributed Space-Time Block Coding", 
IEEE Trans, on Wireless Commun., vol. 54, no. 7, pp. 1195-1206, Jul. 
2006. 

[4] Yulong Zou, Yu-Dong Yao, Baoyu Zheng, "Opportunistic Distributed 
Space-Time Coding for Decode-and-Forward Cooperation Systems", 
IEEE Transactions on Signal Processing, vol. 60, pp. 1766 - 1781, 
April 2012. 

[5] J. Abouei, H. Bagheri, A. Khandani, "An Efficient Adaptive Distributed 
Space-Time Coding Scheme for Cooperative Relaying", IEEE Trans, 
on Wireless Commun., vol. 8, Issue: 10, pp. 4957-4962, October 2009. 

[6] M. Kobayashi, G. Caire, N. Jindal, "How Much Training and Feedback 
are Needed in MIMO Broadcast Channels?", IEEE International Sym- 
posium on Information Theory, 2008. ISIT 2008., p.p. 2663 - 2667, 
6-11 July 2008. 

[7] A.D. Dabbagh, D.J. Love, "Feedback Rate-Capacity Loss Tradeoff for 
Limited Feedback MIMO Systems", IEEE Transactions on Information 
Theory, vol. 52, p.p. 2190-2202, May 2006. 

[8] J. Akhtar, D. Gesbert, "Extending Orthogonal Block Codes with Partial 
Feedback", IEEE Transactions on Wireless Communications, vol. 3, 
p.p. 1959-1962, Nov. 2004. 

[9] Choi Ilhwan, Kim Jong-Kyu, Lee Heunchul, Lee Inkyu, "Alamouti- 
Codes Based Four-Antenna Transmission Schemes with Phase Feed- 
back", IEEE Communications Letters, vol. 13, p.p. 749-751, Oct. 2009. 



[10] G. Jongren, M. Skoglund, "Quantized Feedback Information in Orthog- 
onal Space-Time Block Coding", IEEE Transactions on Information 
Theory, vol. 50, p.p. 2473-2486, Oct. 2004. 

[11] DJ. Love, R.W. Heath, Jr., "Limited Feedback Unitary Precoding for 
Orthogonal Space-Time Block Codes", IEEE Transactions on Signal 
Processing, vol. 53, p.p. 64-73, Jan. 2005. 

[12] P. Clarke, R. C. de Lamare, "Joint Transmit Diversity Optimization 
and Relay Selection for Multi-relay Cooperative MIMO Systems Using 
Discrete Stochastic Algorithms", IEEE Communications Letters, vol. 
15, p.p. 1035-1037, Oct. 2011. 

[13] P. Clarke and R. C. de Lamare, "Transmit Diversity and Relay Se- 
lection Algorithms for Multirelay Cooperative MIMO Systems" IEEE 
Transactions on Vehicular Technology, vol.61, no. 3, pp. 1084-1098, 
October 2011. 

[14] R. C. de Lamare, "Joint Iterative Power Allocation and Linear Interfer- 
ence Suppression Algorithms for Cooperative DS-CDMA Networks", 
IET Communications, 2012. 

[15] S.M. Alamouti, "A Simple Transmit Diversity Technique for Wireless 
Communications", IEEE Journal on Selected Areas in Communica- 
tions, vol. 16, p.p. 1451 - 1458, Oct. 1998. 

[16] R. C. de Lamare and A. Alcaim, "Strategies to improve the performance 
of very low bit rate speech coders and application to a 1.2 kb/s codec," 
IEE Proceedings- Vision, image and signal processing, vol. 152, no. 1, 
Feb. 2005. 

[17] R.C. de Lamare, R. Sampaio-Neto, "Minimum Mean-Squared Error 
Iterative Successive Parallel Arbitrated Decision Feedback Detectors 
for DS-CDMA Systems", IEEE Trans, on Commun., vol. 56, p.p. 778- 
789, May 2008. 

[18] S. Haykin, "Adaptive Filter Theory", 4th ed. Englewood Cliffs, NJ: 
Prentice- Hall, 2002. 

[19] R.C. de Lamare, R. Sampaio-Neto, "Blind Adaptive MIMO Receivers 
for Space-Time Block-Coded DS-CDMA Systems in Multipath Chan- 
nels Using the Constant Modulus Criterion", IEEE Trans, on Commun., 
vol. 58, no. 1, Jan. 2010. 

[20] B. Sirkeci-Mergen, A. Scaglione, "Randomized Space-Time Coding 
for Distributed Cooperative Communication", IEEE Transactions on 
Signal Processing, vol. 55, no. 10, Oct. 2007. 



